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Experimental and theoretical study of superconductor (S) - ferromagnet (F) - superconductor junctions 
showed that in certain range of parameters (e.g., the length of the ferromagnet dp, the exchange field, -Bex) the 
ground state of a SFS junction corresponds to superconducting phase difference n or 0. The phase diagram 
of a SFS junction with the respect to n and states is investigated in this letter in Eex,dF,T space. It is 
shown that the phase diagram is very sensitive to the geometry of the system, in particular, to the amount of 
disorder disorder. 
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PACS: 74.50.+r, 74.80.-g, 75.70.-i 

Recently many interesting phenomena were investi- 
gated in Superconductor (S) - Ferromagnet (F) - Super- 
conductor Josephson contacts. One of the most interest- 
ing effects is the so called 7r-state of SFS junctions [1-11] 
in which the equihbrium ground state is characterized 
by an intrinsic phase difference tt between the two su- 
perconductors. 

Theoretical study of SFS junctions [1-3] showed that 
if Ecx is fixed in the ferromagnet then the 7r-state ap- 
pears at < dp^ < dp < d'^^ 
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the critical current 



etc., and near dp^ (i 
Icidp) has a cusp. Recent experiments [4, 5] showed 
that the critical current - temperature curve, Ic{T), in 
SFS junctions at dp ~ dp\ i = 1,2 has also a cusp. The 
temperature of the cusp was identified with tt — tran- 
sition temperature. These experiments became the mo- 
tivation of the theoretical investigations of Ic{T) curves 
and phase diagrams of SFS junctions. It was shown in 
Refs. [8, 9] that if dp w d'^p (i = 1, 2, . . .) in very short 
ballistic SFS junctions {dp < Co = hvp/A{T = 0)) then 
there is tt — transition at certain temperature and the 
TT state is always (for all i) at larger temperatures then 
0-state. This prediction is in contradiction with the ex- 
perimental data and calculations of the phase diagram 
in dirty SFS junctions based on linearized Usadel equa- 
tions, see Ref. [5] and Refs. there in, where the order 
of TT and phases with the respect to the temperature 



d'-pK SFS 



is one at dj? « d^p^ and the opposite at d^ 
junctions investigated in Ref. [5] were dirty. From the 
first glance it may seem that disorder strongly influences 
on the phase diagram of SFS junctions. This is exactly 
so. It is shown in this letter that the phase diagram 
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of SFS junctions is rather sensitive the geometry of the 
system, in particular, to the amount of disorder in the 
junction. 

The paper is organized as follows: most of the paper 
is devoted to investigations of phase diagrams of SFS 
junctions near the first cusp of Ic{dp) (at dp « d^-*), 
and the case d_F « dp\ i > 1 is discussed in end. 

The superconducting SFS junction that is investi- 
gated here is sketched in FigQ A barrier (e.g., an in- 
sulator layer) is situated at position x — a from the 
junction center. 

First will be considered ballistic SFS junctions. I 
assume that the exchange energy of the ferromagnet 
-Eex Ep; there is no barrier at SF boundaries: the 
probability of Andreev reflection of subgapped Bogoli- 
ubov quasiparticles from a SF boundary is equal to 
unity; variation of the superconducting gap in S near 
the boundary will be neglected (this is correct approxi- 
mation if our SFS is a quantum point contact [12,13]; in 
other cases this approximation can be used because it 
usually leads only to few percent mistakes in the Joseph- 
son current). Then the Josephson current can be found, 
for instance, using scattering matrix method [8, 12]: 
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where 



g{Lu, ip, a) = {2lo^ + A^) cosh($) + 



-I- 2tj\/a;2 -)- A2 sinh($) -f T cos ip + TZ cosh (3 . (2) 

Here ip is the phase difference between the superconduc- 
tors; A is the bulk superconducting gap; w = 27r(n -|- 
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Fig.l. A sketch of a SFS junction. A barrier (e.g., an 
insulator layer) is situated at position x = a from the 
junction center. 



1/2), n = 0,±1, . . .; $ = 2dF{i^ + iEc^(T)/hvF cos e, {3 = 
Aa(uj + iEcyia)/hvF cosB and T = TZ = r^ri, where 
(i|)^, (rj)^ are transmission and reflection probabilities 
of the barrier for spin- up electrons, and 9 is the angle 
between the trajectory and the AT-axis. Eqs.QJ-© can 
be generalized if the ratio E^-^/Ep is arbitrary. Then, 
for example, $ = lmdF{\/2m{EF + iuj + E^^a) — 
yj2m{EF — iuj — ii^exc)}, P can be written in a simi- 
lar way. w-dependence of <& was neglected in Refs. [7-9] 
because there dp was much smaller then ^o- Physically 
ballistic model of a SFS junction could be reahzed, for 
example, in gated heterostructures [14] or in the break- 
junctions [15] in external magnetic field producing Zee- 
man splitting of Andreev levels. [16] 

The case a = (then the scattering potential of the 
junction is symmetric), d_F ^ C = ^^f/A was considered 
in the papers [7-9] . This model is important because on 
qualitative level it well describes SFS junctions where F 
is a ferromagnetic granule or a spin-active interface, see 
Ref. [7] and refs. therein. It was shown that if one fixes 
Eax and changes the temperature then the 7r-state of a 
ballistic SFS junction appears usually at higher temper- 
atures then 0-state. This is not so if the junction is dirty 
and there is no other scattering potential in the junction 
then the impurity potential. To illustrate this point of 
view i consider Josephson current in two similar SFS 
junctions with the same dimensionless normal conduc- 
tance per channel, same A, i?ox, number of open chan- 
nels, etc... The only difference between the junctions 
that in the first one the normal conductance is provided 
by an insulator layer with a flat surface at the center of 
the ferromagnet like in Fig^and in the second one — 
by nonmagnetic impurities in the ferromagnetic region. 
Following parameters were chosen: dp = £,o = hvp/'^o, 
Ao = A(T = 0); e = 2dp'i;cx/7i-CoAo = 0.18. The dis- 
order strength in the second junction was Aq/t = 10, 
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Fig. 2. The figures show critical current - tempera- 
ture relations in two short SFS junctions with the same 
normal conductance, near the first minimum of Iddp) 
(at dp ~ d'^p''). Figl^Ji corresponds to ballistic junc- 
tion and Fig|2l3 — to dirty one. The insets show 
"phase diagrams" of the junctions. The current scale 
Jo — A'^chcAo/ft, where A'^ch is the number of open chan- 
nels in the junction. The thick curve in Fig|5Ji corre- 
sponds to the same O — 0.18 as in Fig^; the other 
curve in Fig|5^ corresponds to O « 0.5. It is seen that 
the TT-state in a ballistic SFS junction is at higher tem- 
peratures then 0-state; the opposite phenomenon takes 
place in a dirty SFS junction. 



where rvp is the mean free path. The transmission 
probability of the insulator layer in the first junction was 

D{e) = i:)ocos2 6'/(i-i:)o-i-i:'ocos2 6'), with Dq = 0.127. 

The current scale Iq = NcheAo/h, where Net is the 
number of open channels in the junction. The thick 
curve in Fig|2^ corresponds to the same Q = 0.18 as 
in Fig^; the other curve in Fig|2K corresponds to 
« 0.5. It follows that near the first cusp of Iddp) 
(i.e., dp ~ d^F^) TT-state in a short ballistic SFS junction 
is at higher temperatury then 0-state; the opposite phe- 
nomenon takes place in a dirty SFS junction. If there 
is an insulator barrier in a ballistic short SFS junction 
and we add nonmagnetic impurities in it then the tran- 
sition from the case shown in Fig|2Ji to the case shown 
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Fig. 3. The figure exhibits critical current - temperature 
relation in a long d/ ~ 3^o ballistic SFS junction with 
Do = 0.127, Q ~ 0.445, near the first minimum of 
Icidp) (at dp ~ dp^). The inset shows the "phase 
diagram" of the junction. If impurities were added in 
this type of junctions then the phase diagram would fi- 
nally look like in Figl^J). If i make Do ~ 1, see Fig0] 
then the phase diagram will become similar to the phase 
diagram shown in Figl^Jj. 



in Figl^b is expected to occur when I /dp ^ D, where 
I is a mean free path. Josephson current calculations 
in ballistic junctions were performed using Eqs.(^-lI2Jl 
with t| = t| = in dirty junctions — using Ricatti 

form of Eilenberger quasiclassical equations [17]. 

A "short" SFS junction with dp ^ was consid- 
ered above. What may happen if dp > is illustrated 
in Fig|3| The junction is ballistic (with the same nor- 
mal conductance corresponding to Dq — 0.127) and 
dp « 3^0- The inset shows the "phase diagram" of 
the junction. So by freezing long SFS junction (at 
dp « dp"^) starting from Tc one can go through a se- 
quence of phase transitions: ^ tt ^ 0. In general one 
can make junction phase diagram similar to any phase 
diagram depicted in Figs lll3l by the proper choice of 
the position of the insulator barrier and the length of 
the junction. For example, if i make Dq — 1 then the 
phase diagram in Fig^will transform to the phase dia- 
gram shown in Fig|2b (the junction is ballistic!). This is 
shown in Fig01 The same transformation of the phase 
diagram occurs due to impurities as show numerical cal- 
culations. 

Below is given a short description of numerical cal- 
culation procedure that was used for drawing of Fig|2|3. 
Calculations were done using Ricatti representation of 
the Eilenberger quasiclassical equations (because Eilen- 
berger equations are unstable) [17]. The quasiclassical 
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Fig. 4. The figure shows critical current - temperature 
relation in a long d/ ~ 3^o ballistic SFS junction with 
Do = 1, O « 0.42, near the first minimum of Iddp) 
(at dp ~ dp^). The inset shows the "phase diagram" of 
the junction. The phase diagram (and the current tem- 
perature relation) in a short dirty SFS junction, Fig|5|3, 
looks quite similar to one shown in Fig|S] 



Green functions can be parameterized via the new func- 
tions a and b, so 

2a . 1-ab . 

^^TT^^^s""' ^^TT^^^sn-. (3) 

The amplitudes a, b change according to the Ricatti 
equations 

vpVa + 2ujjia + Kro^ ~ Ar = 0, (4) 
vpVb - 2LURb - Anb'^ + Ar = 0, 

where Ar ^ A + Ar ^ A* + ^ and ujr = 
uJn + ifjEc^ + {g)/2T. Here (...) denotes averaging 
over directions of the quasiclassical trajectories, r is the 
mean free path of an electron in the impurity poten- 
tial. Eqs.Q) were solved numerically self consistently, 
the Josephson current density was evaluated as follows: 

j = -ineiyT^^{vpg)^^, (5) 

ui a 

where v is the normal density of states. I tried to ex- 
pand the Ricatti equations over t^^ in the first order to 
find analytically how disorder influences on the phase 
diagram of a short SFS junction. Calculations showed 
no effect in the first order. 

The discussion above may lead to a conclusion that 
in dirty Josephson junctions with ferromagnet layers be- 
tween superconductors the phase diagram always looks 
similar like in FigI23 if dp ~ dt)^ . It seems that this 



4 



N.M. Chtchelkatchev 



statement is true for short SFS junctions, but this is 
not true in general. To prove it i shall give below an 
example (more examples and detailed discussion will be 
give in the extended version of this letter). Consider a 
dirty SFIFS junction (e.g., like in FigQJi). The Joseph- 
son current can be found from Usadel equations [17]. In 
general Usadel equations are nonhnear over quasiclassi- 
cal greens functions. But near the critical temperature 
of the junction or if superconductors and ferromagnets 
are weakly coupled the Usadel equations can be lin- 
earized. Then the linearized Usadel equations in the 
ferromagnets for anomalous function / look like: 



2\lu\ + 2i Ji. 2(7 sign 

1^2 



/ = o. 



(6) 



where Di^2 is the diffusion constant and Ji_2 — the ex- 
change energy in the left (right) ferromagnetic layer. 
If the magnetizations of the ferromagnetic layers are 
collinear then the /-function to the left of the insula- 



tor layer I, see Fig^, is / = 
and to the left — / 



\/\i^iIpi\ 



-(Ce''^^ 



(^e«i^ + Be"«i^) 
De-''^^). Here 



Ki,2 = \/(|w| -I- HJx^2<y sign a;)/Di,2, where p\,2 is the 
resistivity of the ferromagnets. The amplitudes A, B, . . . 
are not independent but they are connected by bound- 
ary conditions [18] that I write here in the matrix form: 
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Fig. 5. The figure shows critical current - temperature 
dependence in the dirty SFIF'S junction with parame- 
ters: Ji = -l.QTc, J2 = 3Tc, dp/^ = 1, a/C = -0.2, 
Di = D2, (^ 



^/Djn), R/iPiO = 0.1. Here dp is 
The critical current is normalized on 
its value in the maximum of Ic(T) corresponding to 0- 
state. 



agam near dp 



Aexp(±i0/2) are the gaps of the left 



R 



where Al(7?) 

(right) superconductors, Q. 

the resistances of the SF boundaries. From Eas.(f7|l- Hll|) 
i get: 



Mr) 



s = 



K,i-^K,2 — K,iK.2R 



^(k.i — K.2)a V K,lK2 

k2-\-P^i — kik2R 




(12) 



where S is the "scattering matrix" of the F-F bound- 
ary (diagonal elements of S play the role of "reflection" 
amplitudes, off diagonal — "transmission" amplitudes). 
Here ki^2 — 1^1,2/91,2] R is resistance of the F-F bound- 
ary. When R — the scattering matrix S is similar to 
the quantum mechanical scattering matrix of a poten- 
tial step. At the SF interface the boundary conditions 
look like: 



(9) 
(10) 
(11) 




The approach applied above is similar with the scat- 
tering matrix method used in Ref. [12] for calcula- 
tion of the Josephson current in SNS junctions. Then 
the Josephson current density can be found: j — 
^r^T.^ ^^{f9f - fdf}, where Tr is taken over spin 
degrees of freedom and /(w) — f*{—uj). Fig [S] shows 
critical current - temperature dependence in the dirty 
SFIF'S junction with parameters: Ji = — l.QTc, J2 — 
STe, dp/^ = 1, a/^ - -0.2, - D2, (e = yWTc), 
R/iPiO = O-l- Linearization of the Usadel equations 
is correct at these parameters. Figgis similar to FigEl 
corresponding to ballistic SFS junction. 

The phase diagrams were considered above in 
(T, 6) ~ {T,Eex) space at fixed dp. If one fixes i?ex 
near the first cusp of Ic{Ecx) and changes dp he will 
obtain similar figures. 

Everywhere above I considered phase diagrams of 
SFS junctions near the first cusp of Ic{dp) (at dp w 
d^'p). Below I briefly discuss the phase diagram in gen- 
eral case. FigEl shows tree phase diagrams. In all the 
figures 2Ecx/t^Ao — 1- The first and the second di- 
agrams, FigsEIa-b), correspond to ballistic SFS junc- 



Transitions between "tt" and "0" states in superconductor - ferromagnet - superconductor junctions. 5 



1.0 ■ 
0.8 
0.6 
0.4 
0.2 
0.0 
, 0.8 
0.6 
0.4 
0.2 
0.0 
0.8 
0.6 
0.4 
0.2 



- a) 

J 0-state 
, 1 , 1 


/ 

' K-state / 
1,1,1,1,/, 


D-state 


- b) 

" 0-state 


\ Ti-state / 
1 , N \ . 


0-state 


1 0-state 
' 1 ' 


7i-state / 0-state 
1 1 1 i 'i i 1 i 1 i 1 i 



0.0 0.2 0.4 0.6 0.8 



1.0 1.2 1.4 1.6 



F ^0 

Fig. 6. This figure illustrates how geometry of a SFS 
junction and disorder could influence on the shape of 
the phase diagram. In all figures 2i?ex/7rAo = 1. The 
first and the second diagrams, Figs|SIa-b), correspond 
to ballistic SFS junctions described by Eqs.(0-l|5J with 
Do = 1 (no layer I, see Fig^J and Do = 0.127, a = 
correspondingly. The role of impurities can be seen 
in Figinic). There Do = 1 and h/{AoT) = 10. At 
dp/^o ~ 1 the normal conductance of the junction in 
FiglSJc) becomes equal to the normal conductance of 
the junction shown in Fig|^a). If I draw the phase di- 
agrams for the SFS junctions in Q,Tc space it will look 
similar to Fig|S] 



tions described by Eqs.l^-lO with Dq = 1 (no layer I, 
see Fig^l and Dq — 0.127 correspondingly. The role 
of impurities can be seen in FigEfc). There Do = 1 
and h/{AoT) = 10. At dp/£^Q = 1 the normal conduc- 
tance of the junction in Fig|HIc) becomes equal to the 
normal conductance of the junction shown in Fig|SJa). 
If I'd use the linearized Usadel equations to describe 
the phase diagram of a SFS junction then I'd get the 
graph like in FigE^c). The phase diagram in Fig|SJc) 
qualitatively agrees with experimental results in Ref. [5] 
that the order of tt and phases with the respect to the 
temperature is one at dp w d^p^ and the opposite at 



dp 



7(2) 



In conclusion, tt — transitions in Josephson junc- 
tions with ferromagnetic layers are investigated in this 
letter. It is shown that the phase diagram is very sensi- 
tive to the geometry of the system, in particular, to the 
amount of disorder in the junction. 
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